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ABSTRACT 



O We consider the recent limits on dark matter - nucleon elastic scattering cross section from the 
Q ■ analysis of CDMS II collaboration using the two signal events observed in CDMS experiment. 
•• . With these limits we try to interpret the Super-Kamiokande (SK) bounds on the detection rates 
• i-h ', of up-going muons induced by the neutrinos that are produced in the sun from the decay of 
^ | annihilation products of dark matter (WIMPs) captured in the solar core. Calculated rates of 
d ' up-going muons for different annihilation channels at SK using CDMS bounds are found to be 
orders below the predicted upper limits of such up-going muon rates at SK. Thus there exists room 
for enhancement (boost) of the calculated rates using CDMS limits for interpreting SK bounds. 
Such a feature is expected to represent the PAMELA data with the current CDMS limits. We also 
show the dependence of such a possible enhancement factor (boost) on WIMP mass for different 
WIMP annihilation channels. 
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Weakly interacting massive particles (WIMP) as dark matter in our galaxy can be trapped 
inside massive heavenly bodies like sun due to later's gravity [1]. This gravitational trapping may 
happen when such dark matter in course of its passage through sun undergoes elastic scattering 
off the nuclei present in the solar core which causes its final velocity to fall below its velocity 
of escape from solar gravitational pull. These trapped dark matter may undergo the process of 
pair-annihilation producing primarily b, c and t quarks, r leptons, gauge bosons etc. Mass and 
composition of dark matter determine the annihilation products which in turn produce neutrinos 
and antineutrinos either through decay or pair annihilation. Such neutrinos from the dark matter 
annihilation in solar environment have also been studied by previous authors (e.g. [2]). Detection 
of such solar neutrinos in terrestrial detectors not only provides indirect evidence of dark matter 
but along with the results from the direct detection experiments of dark matter such as CDMS 
[3], DAMA [4], XENON [5] etc. provides more insight into the nature of dark matter and its 
interactions. Analysis of recent observations of two dark matter signal events with the Cryogenic 
Dark Matter Search experiment (CDMS II) at the Soudan Underground Laboratory combined 
with all previous CDMS II data has been reported by CDMS collaboration [6]. It sets new 
upper limits on the WIMP-nucleon elastic scattering cross section (a x ) as a function of WIMP 
mass (m x ) [6]. The indirect searches for WIMPs through their annihilation in sun, with 1679.6 
live days of data from SK detector using neutrino-induced upward through-going muons provide 
WIMP-induced upward muon flux limits at SK as a function of WIMP mass [7]. In this work we 
use the 90% Confidence Level (C.L.) limits on a x (m x ) from recent CDMS analysis [6] to calculate 
corresponding limits on detection rates of up-going muons at SK as a function of WIMP mass 
(m x ) and compare them with the results in [7]. 

The differential flux of neutrinos of type i (i — z/ M , v^) at earth from WIMP annihilation 
products in the sun is given by [8] 
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where R is sun-earth distance and Bp is the annihilation branch for channel F. (dN/dE) Fi is the 
differential spectrum of neutrinos of type % in the sun for the annihilation channel F. The total 
rate for WIMP annihilation in the sun, Ya is given in [10] 

T A = ^tanh 2 [(aC) 1/2 r] (2) 

where r ~ 4.5 Gyr is the age of sun. a = {crv)/A\/2V is a function of the average WIMP 
annihilation cross section and the effective volume V of WIMPs in the sun [10, 11, 1, 12, 13, 14]. 
Under the astrophysical assumptions on density and velocity distribution as mentioned in [10, 11] 
(local dark matter density, p x = 0.3 GeV cm~ 3 , mean velocity of dark matter, v = 300 km 
Maxwellian distribution of velocities etc.) the dark matter capture rate C in the sun is 
approximated as a function of the ratio of the WIMP-nucleus elastic scattering cross section to 
the dark matter mass as [10, 1] 

C ~ 10 29 ^ GeV pb^sec- 1 (3) 
m x 
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For such indirect detections of WIMPs at SK, neutrinos are detected through up-going muons 
produced by charged current interactions of neutrinos with the rock below the detector. With 
numerical values of this cross section, the muon range in the rock and expressing the energy 
distribution of neutrino flux in terms of its second moments, the total detection rate at SK detector 
of up-going muons induced by neutrinos from WIMP annihilation in the sun is given by [8] 



detect 



(1.27 x lO^yr- 1 )^ (£v) 2 E «.^E^(iVzV, x A eff (4) 



i=v,u F 

where, A e ff(m 1200m 2 ) is the muon effective area of the SK detector [15], a^s are the neutrino 
scattering coefficients. The range of neutrino induced muons in the rock are given as the coefficients 
These coefficients are given by a v = 6.8, a v = 3.1, b u = 0.51, bp = 0.67 [8]. (N z 2 ) Fi ^ are 
the second moments of the spectrum of neutrino type i for the WIMP annihilation channel F in 
the sun. The {N z 2 ) Fi for different channels relevant for the present calculations are listed below 
[8,9]. 

(a) rf channel: 

(Nz 2 )^^) = T T ^ uD h T)i (E- m jT t ) (i = u II ,u fl ) (5) 

where E m j is the injection energy of the decaying WIMP annihilation product inside the sun, 
the branching ratio T^^p ~ 0.18, t u (t d ) = 1.01 x 10~ 3 (3.8 x 10~ 4 ) GeV -1 are the stopping 
coefficients for v{y) and 

4 + 2/ 
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(b) bb channel: 



168 + 354y + 348y 2 + 190y 3 + 56y 4 + 7y 5 
K ^ V) = 1260(1 + y) (6) 



(N z \(E inj ) ~ T^x^h^^Ef)^ (7) 



where, the branching ratio T^^x = 0.103. The hadronization and the decay processes of 
the quarks from WIMP annihilation in the sun are characterized by the mean energy (Ed) 
of the hadron, (E d ) = E c exp (f|) E 1 (f|), E c = 470 GeV [9] . E = Z f E inj is the initial 
hadron energy for quarks injected with energy E in j and Zf(— 0.73) is the quenching fraction 
for 6-quarks to account for the loss of energy during hadronization, E 1 (x) = ^j-dy. Also 
(Efj = E C (E - (E d )) and h b)i is same as h T>i . 



(c) W + W~ and ZZ channel: 

(Nz 2 ) t (E m j) 



" w ^ u 2 + 2ET t (l + a t ) + E 2 ria l (l + a t " 1 ' '"■ ' " ' 2 



2T Z _ U ^ 2 + 2Er l (l + a t ) + £ 2 rf a,(l + a,) 



E=E ini (l+P)/2 
E=E inj (l-P)/2 



fi El-rfa^ - 1)(1 + En)-^ 



(8) 
(9) 



£=£ inj (l+/3)/2 
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Figure 1: CDMS upper bounds on WIMP induced muon event rates in different annihilation 
channels with E in j = m x (left panel) and = m x /3 (right panel) 

where, a u (a p ) = 5.1(9.0) and the branching ratios Tw^v — 0.105 , Tz-+vupu = 0.067. (3 
being the velocity of the gauge bosons. 

The injection energy E- 1X ^ is the energy with which the WIMP annihilation products bb, tt, W 
and Z etc. are produced. In this work we present our results for two benchmark scenarios namely 
-Einj = Tn x an d 1^- 

Recently the CDMS collaboration announced two dark matter signal events with 90% C.L.[6] 
which set new limits on dark matter-nucleon scattering cross-section a x for different m x s [6]. In 
this work we choose our set of cr x 's and corresponding m x s from this limit (exclusion plot for 
"Soudan(all)" in Fig. 4 of Ref. [6]). With these CDMS limits on a x we compute the detection 
rates for neutrino induced muons (from the product of the Dark Matter annihilation in the sun) 
for different annihilation channels. Thus we calculate the CDMS induced limits on event rates 
in SK. We compute the detection rate for neutrino induced muons (from the product of WIMP 
annihilation in the sun) using Eq. 4. 

In this work we present our results for two benchmark scenarios namely E m j = m x and A 
conservative estimate for neutrino production from WIMP annihilation in the sun is obtained by 
considering the dominating effect of bb production in the region mf, < m x < mw, that of rf and 
W, Z production in the regions mw < < m t anc ^ m x > m * respectively [8], where m&, m t , mw 
are respective masses of 6-quark, t-quark and ly-boson. The estimated limits of WIMP induced 
muon detection rates at SK using Eqs. (4-9) in the annihilation channels bb ,rf and W, Z as a 
function of WIMP mass has been shown in Fig. (1). These rates are calculated for E in j = m x 
(left panel of Fig. 1) and E- m j = m x /3 (right panel of Fig. 1). Plots for each channel in Fig. (1) 
are obtained by taking the value of branching fraction (Bp) for that channel to be 1. This can 
be justified by the fact that each channel is effective in a particular m x region as discussed above. 
However, we have shown in our results broader m x ranges for each of the channels. Also shown in 
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Figure 2: Plot of B vs m x in each channel with E in j = m x (left panel) and = m x /3 (right 
panel) 

this figure is the corresponding SK limit [7] and our results indicate existence for room for an m x 
dependent enhancement of order 10 3 — 10 5 of the detection rate in order to interpret SK bounds 
in terms of the CDMS a x (m x ) limit. 

In a recent work, Cao et al [16] have shown that to study the anti-proton fraction of the results 
of satellite borne PAMELA experiment [17, 18] using the recent CDMS results one needs to invoke 
a boost factor (~ 10 3 ) for the annihilation of WIMPs into quarks. This, in general is consistent 
with the enhancement needed to explain the excess of positron fraction as observed by PAMELA 
experiment [17]. Our results also indicate similar trends that if one uses the CDMS limit as is 
done in the work, one is allowed to enhance the detection rate for the present of SK bounds. 
We redefine Bp in our framework (Eq. 4) as B x B' F where B' F is the branching fraction for a 
particular process and B is some boost factor. For a single channel to represent the entire neutrino 
production process, Bp = 1. Then we estimate the upper limits of B's for different dark matter 
annihilation channels in the sun for the recent a x — m x CDMS limit. We do this by comparing the 
detection rates of neutrinos - obtained from the decays of the dark matter annihilation products 
in the sun - at Super-Kamiokande (SK) detector. The upper limits on WIMP induced up-going 
muons at SK as a function of WIMP mass is given by the SK Collaboration [7] . In Fig. 2 (left 
panel) we show the variation of B with m x for E^ = m x . Right panel of Fig. 2 shows similar plots 
for E- m j = m x /3. From Fig. 2 (left panel) we see B varies (with m x ) between ~ 10 5 — ~ 10 4 for bb 
channel, ~ 10 4 — ~ 10 3 for rr channel, ~ 10 3 for each of W and Z channels etc. For E in ^ = m x /3 
case (right panel of Fig. 2) however, the variation of -B's for each of W and Z channels is between 
~ 10 3 — ~ 10 2 and for rr channel it is from ~ 10 4 — 10 2 . These limits are consistent with the 
order 10 3 as discussed earlier. 

Although in Figs. (1-2) we have shown the results for W and Z channels separately but as 
discussed earlier both are dominant in the region m x > m t . Therefore we consider the effect of 
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Figure 3: Iso m x plot in B x Bp (W channel) - B x Bp' (Z channel) plane interpreting SK limits 
with E in j = m x (left panel) and E in j = m x /3 (right panel). 

both of them together in calculating the rate using Eq. (4). In Fig. 3 (left panel) we present in 
B x Bp' (W channel) - B x Bp' [Z channel) plane the iso m x plots for E iQ j = m x that represent 
the required rate to interpret the SK limits [7] . Right panel of Fig. 3 gives the similar plots but 
with E [n j = m x /3. 

In the present work we compute the WIMP induced upward going muon rates at SK. They are 
induced by neutrinos from WIMP annihilation products in the sun. We use recent CDMS bounds 
on WIMP-nucleon scattering cross sections for different WIMP masses in our rate calculation. It 
is observed that representation of SK upper bounds on WIMP induced up-going muon rates allows 
an enhancement in the calculated rates in all individual channels. We made an estimation of this 
enhancement as a function of WIMP mass assuming branching fractions for each different channels 
to be 1 (maximum). Interpretation of PAMELA data in terms of CDMS limits also demand such 
enhancements of WIMP annihilation. 
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